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ABSTRACT 

We show that we can obtain a good fit to the present day stellar mass functions (MFs) of a large sample of 
young and old Galactic clusters in the range 0.1 - 10 M© with a tapered power law distribution function with 
an exponential truncation of the form dN/dm oc m" [i -^-('"/'"r) ]. The average value of the power-law index 
a is ~ -2, that of /3 is ~' 2.5, whereas the characteristic mass nic is in the range 0.1 -0.8 Mq and does not 
seem to vary in any systematic way with the present cluster parameters such as metal abundance, total cluster 
mass or central concentration. However, nic shows a remarkable correlation with the dynamical age of the 

3/4 

cluster, namely nic/MQ ~ 0.15 + 0.5 x t^', where Tdyn is the dynamical age taken as the ratio of cluster age 
and dissolution time. The small scatter seen around this correlation is consistent with the uncertainties on the 
estimated value of Tdy„. We attribute the observed trend to the onset of mass segregation via two-body relaxation 
in a tidal environment, causing the preferential loss of low-mass stars from the cluster and hence a drift of the 
characteristic mass iric towards higher values. If dynamical evolution is indeed at the origin of the observed 
trend, it would seem plausible that high-concentration globular clusters, now with median nic — 0.33 Mq, were 
born with a stellar MF very similar to that measured today in the youngest Galactic clusters and with a value 
of wjc ~ 0.15 M0. This hypothesis is consistent with the absence of a turn-over in the MFof the Galactic bulge 
down to the observational limit at ^ 0.2 Mq and, if correct, it would carry the implication that the characteristic 
mass is not set by the thermal Jeans mass of the cloud. 

Subject headings: globular clusters: general — open clusters and associations: general — stars: luminosity 
function, mass function 



1. INTRODUCTION 

There is general consensus that stars do not form in iso- 
lation but rather from the fragmentation of molecular clouds 
that leads to star clusters (e.g. Lada & Lada 2003; Elmegreen 
2010). This makes stellar clusters ideal places to study the 
properties of star formation and its outcome, the stellar ini- 
tial mass function (IMF). However, it is also equally well es- 
tablished today that a large majority of stars, in our Galaxy 
and elsewhere, are not in clusters but in the field, since clus- 
ters disrupt over time (see e.g. Gieles 2010). Therefore, any 
attempt to set constraints on the star formation mechanisms 
from the analysis of the present day stellar mass function 
(MF) of but the youngest clusters cannot ignore the conse- 
quences of their dynamical evolution. 

This issue becomes particularly important if we want to 
compare the results of star formation in physically different 
environments and with various ages. The obvious example is 
addressing differences in the way low-mass stars (< 1 Mq) 
formed in globular clusters (GCs), at z ~ 5 or about 12Gyr 
ago, and in young clusters (YCs) in the local universe. While 
in both cases the raw data show a broad plateau in the mass 
distribution, suggesting a characteristic mass of order a few 
tenths of Mq (e.g. Elmegreen et al. 2008), until the effects 
of dynamical evolution are properly understood, no meaning- 
ful conclusion can be drawn either on the uniformity of the 
star formation process over time or on the possible role of the 
environment. In order to cover a wide range of metallicities 
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and initial densities one is forced to compare stellar systems 
that formed a Hubble time apart, such as the YCs and GCs 
mentioned above. Therefore, the universality of the IMF over 
time and location remains as yet an unresolved issue. 

On the other hand, even though the two-body relaxation 
process that governs the cluster's dynamical evolution and 
that leads to preferential loss of low-mass stars is today rather 
well understood (e.g. Spitzer 1987), it is in general not pos- 
sible to roll back the effects of dynamics and derive the IMF 
from the present day MF, since we cannot trace back the tra- 
jectories of stars that have escaped the cluster It is however 
possible and statistically meaningful to study the evolution 
of the stellar MF on a global scale by looking at the differ- 
ences between clusters at different evolutionary stages. This 
requires a homogeneous sample of high quality observations 
of Galactic (open and globular) clusters, treated in a uniform 
way and with reliable errors. Early in this decade, this type of 
homogeneous study became possible for GCs (Paresce & De 
Marchi 2000), mostly thanks to the Hubble Space Telescope. 
In the meanwhile, high quality data have become available for 
YCs as well, mostly from wide field ground based surveys, 
thereby making this study possible on a global scale. 



2. THE TAPERED POWER-LAW 

Since our goal is to detect and quantify changes in the shape 
of the MF in different environments, in particular in GCs and 
YCs, we need a functional form for the distribution of stellar 
masses that is flexible enough to adapt to the variety of ob- 
served MFs and yet simple enough to be described by a small 
number of parameters over a wide mass range. Kroupa (2002) 
has proposed a segmented, multi-part power-law for this pur- 
pose, but we cannot adopt it here because this approach fixes 
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TABLE 1 

The sample of stellar clusters used in this study 



ID 


Name 


a 


/3 


[Mq] 


m range 

[Mo] 


[Fe/H] 


logMiot 

[Mq] 


c 


logr 
[yr] 


log?dis 

[yr] 


Tdyn 


Ref 


1 


pOph 


-1.9±0.2 


1.9±0.2 


0.17 ±0.03 


0.02- 


-7 


0.08 


3.3 




5.7 


8.7 


0.001 


a 


2 


Orion N. C. 


-1.9±0.2 


2.0 ±0.3 


0.19 ±0.05 


0.02- 


-50 


-0.01 


1.8 


0.3 


6.0 


8.7 


0.002 


b,c 


3 


Taurus 


-2.3 ±0.2 


1.8±0.2 


0.51 ±0.05 


0.03- 


-3.5 


-0.03 


1.7 




6.0 


6.3 


0.5 


d 


4 


IC348 


-1.9±0.1 


3.3±0.2 


0.12 ±0.05 


0.03- 


-2.2 


-0.01 


1.8 


0.3 


6.3 


8.7 


0.004 


e 


5 


o-Ori 


-2.0 ±0.1 


3.2±0.2 


0.14 ±0.02 


0.02- 


-2.4 


-0.01 


2.3 


0.6 


6.5 


8.7 


0.006 


f 


6 


AOri 


-1.9±0.1 


2.3 ±0.2 


0.12 ±0.02 


0.04- 


-2.7 


-0.01 


2.8 


0.6 


6.7 


8.7 


0.01 


g 


7 


Chal 


-1.9±0.2 


2.3±0.3 


0.15 ±0.03 


0.01- 


-3.5 


-0.11 


2.6 




6.7 


8.7 


0.01 


h 


8 


IC2391 


-2.1 ±0.1 


2.4 ±0.3 


0.14 ±0.02 


0.05- 


-0.5 


-0.09 


2.2 


0.6 


7.7 


8.8 


0.06 


i 


9 


Blanco 1 


-1.7±0.1 


1.8±0.2 


0.22 ±0.03 


0.05- 


-1.75 


0.14 


3.5 


0.6 


8.1 


9.1 


0.09 


j 


10 


Pleiades 


-2.2 ±0.1 


2.3 ±0.2 


0.27 ±0.03 


0.04- 


-9 


0.03 


2.9 


0.6 


8.1 


9.3 


0.06 


k 


11 


M35 


-1.7±0.1 


2.4 ±0.2 


0.30 ±0.03 


0.09- 


-1.4 


-0.21 


3.2 


0.7 


8.2 


9.4 


0.08 


i 


12 


Coma Ber 


-1.3 ±0.4 


1.5 ±0.3 


0.45 ±0.10 


0.12- 


-1.1 


-0.08 


2.4 


0.4 


8.6 


8.9 


0.48 


b 


13 


Praesepe 


-2.0 ±0.2 


4.0 ±0.5 


0.20 ±0.02 


0.07- 


-1.5 


0.00 


3.3 


0.7 


8.8 


9.5 


0.18 


1, m 


14 


Hyades 


-2.1 ±0.1 


2.8 ±0.2 


0.45 ± 0.05 


0.07- 


-2.7 


0.14 


2.6 


0.4 


8.8 


9.2 


0.37 


n 


15 


NGC 104 


-2.0 


1.8±0.15 


0.33 ±0.02 


0.11- 


-0.7 


-0.76 


5.84 


2.03 


10.11 


10.93 


0.15 


P 


16 


NGC5139 


-2.0 


2.3 ±0.2 


0.33 ±0.02 


0.13- 


-0.7 


-1.62 


6.17 


1.61 


10.06 


10.75 


0.21 


P 


17 


NGC 5272 


-2.0 


2.4 ±0.2 


0.33 ±0.02 


0.16- 


-0.7 


-1.57 


5.65 


1.84 


10.06 


10.92 


0.14 


P 


18 


NGC 6121 


-2.0 


2.8 ±0.25 


0.35 ±0.02 


0.09- 


-0.5 


-1.20 


4.83 


1.59 


10.10 


10.16 


0.88 


P 


19 


NGC 6254 


-2.0 


1.9±0.2 


0.33 ±0.02 


0.13- 


-0.6 


-1.52 


5.00 


1.40 


10.06 


10.39 


0.47 


P 


20 


NGC 6341 


-2.0 


2.0 ±0.2 


0.30 ±0.02 


0.11- 


-0.7 


-2.28 


5.30 


1.81 


10.12 


10.40 


0.52 


P 


21 


NGC 6397 


-2.0 


2.3 ±0.2 


0.33 ±0.02 


0.09- 


-0.55 


-1.95 


4.65 


2.50 


10.10 


10.27 


0.67 


P 


22 


NGC 6656 


-2.0 


2.3 ±0.3 


0.31 ±0.03 


0.10- 


-0.6 


-1.64 


5.46 


1.31 


10.10 


10.62 


0.31 


P 


23 


NGC 6752 


-2.0 


2.6±0.1 


0.34 ±0.02 


0.10- 


-0.6 


-1.56 


5.15 


2.50 


10.07 


10.50 


0.37 


P 


24 


NGC 6809 


-2.0 


2.3±0.3 


0.30 ±0.02 


0.11- 


-0.6 


-1.81 


5.04 


0.76 


10.09 


10.32 


0.58 


P 


25 


NGC 7078 


-2.0 


2.7±0.3 


0.23 ±0.02 


0.15- 


-0.7 


-2.26 


5.75 


2.50 


10.07 


10.93 


0.14 


P 


26 


NGC 7099 


-2.0 


1.6±0.2 


0.25 ±0.02 


0.15- 


-0.7 


-2.12 


5.00 


2.50 


10.11 


10.39 


0.53 


P 


27 


NGC 2298 


-2.0 


3.0 ±0.2 


0.54 ±0.01 


0.20- 


-0.8 


-1.85 


4.50 


1.28 


10.02 


10.25 


0.59 


q 


28 


NGC 6218 


-2.0 


3.2 ±0.4 


0.62 ±0.03 


0.30- 


-0.8 


-1.48 


4.94 


1.39 


10.10 


10.33 


0.60 


q 


29 


NGC 67 12 


-2.0 


3.2±0.2 


0.80 ±0.03 


0.30- 


-0.8 


-1.01 


4.97 


0.90 


10.02 


10.22 


0.63 


q 


30 


NGC 6838 


-2.0 


2.6 ±0.2 


0.60 ±0.02 


0.30- 


-0.8 


-0.73 


4.23 


1.15 


10.14 


10.21 


0.84 


q 



Note. — Table columns are as follows: ID number; cluster name; best fitting value of a and 1 a uncertainty; best fitting value of /3 and 1 cr uncertainty; 
best fitting value of Wc and la uncertainty; mass range over which the TPL fit is performed; cluster metallicity; total cluster mass; central concentration, or 
the logarithmic ratio of tidal radius and core radius; cluster age; estimated time to dissolution; dynamical time, or the ratio of the cluster age and the time to 
dissolution; bibliographic reference for the MF data, as follows. Fo r YC s: ( a) Luhman & Rieke ( 1999); (b) Kraus & Hillenbrand 12007); (c) Slesnick et aQ 
12004); (d) Salas & Cruz-GonzaTe? FMU); (e) 'Muench et al." FMT); (f) "Caballero 1 2008); (g) Banado v Navascues et al. 1 2004a); (h) Luhman 1 2007); (i) 
iBarrado v Navascues et al. l2004b); (j) Moraux et al. (2007); (k) Moraux et al. (2003); (1) Chappelle et al. (2005); (m) Wilham, Rieke & Stauffei ( 1995); (n) 
iBouvier et al. 1 2008). For dense GCs: (p) Paresce & De Marchi ( 2000). For loose GCs: (q) De Marchi, Paresce & Pulone, ( i2Q0 7). For YCs, the metallicity and 
age are from the bibliography in column Ref. and references therein, whereas total mass and concentration are from Piskunov et al. (2008), although the value 
of Mtoi is from Kiiciik & Akkaya (2010) for Taurus, from Hodapp et al. (2009) for cr Ori and from Dolan & Mathieu (2002) for A Ori. For GCs, the metallicity, 
total mass and concentration are from the 2003 revision of the GC catalogue of Hairis (1996), whereas the age is from Forbes & Bridges (2010). The values of 
?dis have been determined as explained in Section 5, whereas r^y^ is the dynamical age of each cluster, expressed as f/fjis- 



the mass points between which the slope is fitted, thereby 
making it impractical or even impossible to detect variations 
in the characteristic mass, which, as we will see, are the sig- 
nature of dynamical evolution. 

Our choice falls instead on a tapered power-law (TPL) dis- 
tribution of the type 



/ (m) = ex m 

dm 



1- 



-(m/mc) 



(1) 



where nic is the characteristic mass, a the index of the power- 
law portion for high masses and /3 the tapering exponent that 
defines the shape of the MF below the characteristic mass nic- 
The TPL combines in one expression the ubiquitous power- 
law shape above 1 M© (e.g. Salpeter 1955) with the plateau 
and drop observed near the H-burning limit (e.g. Paresce & 
De Marchi 2000; Chabrier 2003). As we already showed (De 
Marchi, Paresce & Portegies Zwart 2005), the TPL fits re- 
markably well the MFs of both YCs and GCs over the entire 
stellar mass range. 



3. THE DATA SAMPLE 

The data used in this work come from a number of differ- 
ent sources. As regards GCs, we considered the entire HST 
sample of Paresce & De Marchi (2000), which is made up 
of twelve relatively dense clusters, and added to it four low- 
concentration GCs studied by our team in recent years (for 
details see Table[Tli. De Marchi et al. (2007) showed that a 
strong trend exists between a cluster's central concentration c 
(defined as the logarithmic ratio of the tidal radius and core 
radius) and the shape of its present stellar MF: dense GCs, i.e. 
those with c > L5, have a relatively steep MF with an aver- 
age power-law index of -L5 in the range 0.3- 0.8 M© , while 
loose GCs with c < L5 show a much flatter or even inverted 
MF over the same mass range. Several interpretations exist 
for this result (e.g. Baumgardt, De Marchi & Kroupa 2008; 
Kjuijssen & Portegies Zwart 2009; Vesperini, McMillan & 
Portegies Zwart 2009; Kruijssen 2009), but they all point to a 
shorter time to dissolution for loose GCs. For this reason, it 
is important to cover both dense (c > 1 .5) and loose (c < 1 .5) 
clusters when studying the temporal evolution of the stellar 
MF 
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Fig. 1. — Mass functions of selected clusters in our sample (see Table[T] 
for the ID numbers). The data for YCs are shown as large circles (except 
for Praesepe, where also small circles are used to show the MF at the high- 
mass end as measured by Wilhams et al. (1995); NGC6397 is shown as 
representative for the class of dense GCs (squares); while the MFs of loose 
GCs are indicated by triangles (only NGC2298 and NGC6712 are shown). 
The typical 1 cr uncertainty on the MF is 0.5 dex at ~ 0.3 Mq . The dashed 
lines represent the best fitting TPL distributions, as per the parameters of 
Table[T] with the values of nic indicated by the dark down-pointing triangles. 
For reference, we also show the MF of the Galactic disc. At the low-mass 
end the data are those recently obtained by Covey et al. (2008) from Sloan 
Digital Sky Survey and Two Micron All Sky Survey data (large diamonds), 
whereas above ~ 1 Mq we show the MF obtained by Reid et al. (2002) and 
corrected for the effects of stellar evolution (small diamonds). 



As regards YCs, the data cover high quality observations 
of fourteen star-forming regions and associations of various 
ages, as shown in Table[T] For each YC, we searched for the 
most accurate determination of the MF in the recent literature 
(see references in TablelT]!, looking specifically for informa- 
tion on the number of stars measured in each mass bin. In 
those cases in which the MFs were not directly available in 
tabular form, we extracted the data points (i.e. number of stars 
in each mass bin) from the published graphs, but we specif- 
ically ignored any fits (with power-law or log-normal func- 
tions) or interpolations of these data that the various authors 
may have carried out. The MFs in our sample were deter- 
mined either by converting an observed luminosity function 
via a theoretical mass-luminosity relationship or by counting 
the number of stars falling between evolutionary tracks in a 
color-magnitude or Hertzsprung-Russell diagram. We only 
restricted our search to works in which the MF is based on 
known cluster members and a correction for photometric in- 
completeness has been applied to the luminosity function, if 



needed. Furthermore, to guarantee that the sample is as homo- 
geneous as possible, we have ignored any corrections to the 
MF meant to account for the presence of binaries, since these 
are rather uncertain and only available in a limited number of 
cases. Therefore, the MFs that we consider in this study, for 
both YCs and GCs, are the MFs of stellar systems, i.e. we do 
not distinguish between single and multiple stars. 

In order to prevent biases in the results, the data should re- 
fer to the global MF (GMF), i.e. to the MF of the cluster as 
a whole. For GCs, where complete cluster coverage is not al- 
ways possible, we used information on mass stratification and 
mass segregation to derive the GMF from the local MF (see 
e.g. De Marchi et al. 2006). Alternatively, we used the MF 
measured near the half-light radius (i.e. the effective radius 
containing half of the cluster's luminosity, see e.g. Portegies 
Zwart, McMillan & Gieles 2010), since it has been shown to 
reflect quite reliably the properties of the GMF (De Marchi et 
al. 2000). For YCs, we specifically selected those objects and 
studies for which the coverage is as complete as possible. 

We then performed a multivariate fit to each MF with a TPL 
distribution and derived the values of the parameters a, /3 and 
nic that simultaneously provide the smallest residuals. The 
values of the best fitting parameters and their associated 1 a 
uncertainties are given in Table[T] for all the clusters in our 
sample. The MFs of a few selected clusters are also shown 
graphically in Figure[Tl together with the best TPL fits as per 
the parameters of Table[T] We compare and discuss the results 
of our best fits in Sections, but we precede that with a short 
discussion on the value of a. 

4. MASS FUNCTION SHAPE ABOVE - 0.8 M© 

In the TPL distribution described by Equation 1, the value 
of a determines the shape of the MF above the characteristic 
mass nic- Since the data for GCs do not constrain the mass 
range above the main-sequence turn-off (> 0.8 M0), we had 
to assume a value of a for those objects. An obvious first 
choice would be the slope of the Salpeter (1955) IMF, namely 
a = -2.35. According to Kroupa (2002), that is the average 
value of the slope of the stellar MF in clusters and associations 
above 0.5 M© . However, this is not a meaningful average 
because of the very large scatter at m > 1 Mq displayed by the 
MF slopes included in Kroupa' s (2002) compilation, which in 
turn are taken from Scalo (1998). In fact, after discussing 
with this latter author (J. Scalo, priv. comm.), we realised that 
the a = -2.35 value obtained in that way is not a meaningful 
average. It is also not consistent with any of the individual 
a values that we measure for the YCs in Table[T] except for 
Taurus and possibly Pleiades. 

As an alternative we could use the MF of the Galactic disc, 
but decided against it because of the large uncertainties in- 
volved. The most recent determination is that by Covey et al. 
(2008) using Sloan Digital Sky Survey and 2 MASS data, al- 
though these observations do not reach above ^ 1 Mq (see 
large diamonds in Figure[T]i. To extend the MF to higher 
masses, these authors use the measurements of Reid, Gizis 
& Hawley (2002). We do the same and show with small dia- 
monds in Figure[T]the MF that Reid et al. (2002) obtained us- 
ing an empirical mass-magnitude relationship and after cor- 
rection for the effects of stellar evolution (this provides them 
with an estimate of the IMF). The best fit with a TPL distribu- 
tion gives a = -2.4 ± 0.4, /3 = 2.3 ± 0.4 and m^ = 0.4 ± 0.05, 
revealing again a rather large uncertainty around the value of 
a. 

In the end, following the suggestion of an anonymous ref- 
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eree, we concluded that the most appropriate value of a 
to adopt for GCs is the average of the best fitting a val- 
ues that we obtain in the same mass range for YCs, namely 
a = -1.97 ±0.17. This is also conveniently very close to the 
value that Van Buren (1985) found for the MF slope in the 
range 5- I2OM0 from a detailed analysis of the all-sky sam- 
ple of O stars of Garmany, Conti & Chiosi (1982) that makes 
use of high-resolution extinction maps derived spectroscopi- 
cally to account for the effects of obscuration on star counts 
(Scalo 1986). 

Finally, we note that observational constraints on the initial 
value of a for GCs could in principle come from the anal- 
ysis of the white dwarf luminosity function. The results of 
Richer et al. (2004) on the white dwarf cooling sequence 
of M4 (NGC6121) seem to suggest a value in the range 
-2.5 < a < -1.5, but more recent work by the same team 
(Richer et al. 2008) on another cluster, NGC 6397, seems to 
require a much steeper slope, possibly as large as a ~ -5. 
Therefore, it appears that the uncertainties are presently still 
too large for this method to place meaningful empirical con- 
straints on the value of a for GCs, although progress in this 
field is likely with the advent of the James Webb Space Tele- 
scope in the near future. 

5. EVOLUTION OF THE MASS FUNCTION 

The mean a and (3 values for the YCs in our samplfl 
resulting from our best fits shown in TablefT] are respec- 
tively -1.97±0.17 and 2.51 ±0.62 (la standard deviation 
uncertainties). For GCs, having assumed a = -2, we find 
/3 = 2.49 ±0.47. In spite of the widely different physical prop- 
erties of the clusters in our sample, e.g. in terms of chemical 
composition, total mass, stellar density and age (see Table[Tll, 
the values of a and f3 span a relatively narrow range and, in 
particular, there is no significant difference in /? between YCs 
and GCs. Obviously, having assumed a = -2 for GCs has in 
principle a systematic effect on /3, but this does not seem to be 
large (an even smaller effect on m^. is discussed below). For 
example, with our original assumption of a = -2.35 we found 
(3 = 2.6 ±0.3. This suggests that that there does not seem to 
be substantial difference in the way the stellar MF rises to 
and falls from its peak, i.e. the characteristic mass, in clusters 
covering a wide range of physical properties, including age, 
chemical composition, total mass and density. 

Where large differences are seen, however, is in the value 
of the characteristic mass itself. The average value of nic for 
dense GCs is 0.31 ±0.03 Mq and for loose globulars it grows 
to 0.64 ± O.IOM©, although the statistics here is limited. As 
noted above, the actual value of mc depends in principle on the 
assumed a, but changing the latter by 0.4 would cause m^ to 
vary by just 0.01 M© for the GCs in our sample. It is therefore 
unlikely that the differences in the nic values that we observe 
are just the result of uncertainties in a. The variation in m^ 
is equally large for YCs, ranging from O.I2M0 to 0.51 M0, 
and here the data provide much more solid constraints on the 
value of a. 

In principle, all these m^ variations could reflect differences 
in the star formation process, in turn resulting from differ- 
ences in the physical conditions of the environment (e.g. Lar- 
son 2005 and references therein). However, as the data in Ta- 
ble[T|show, there does not seem to be any correlation between 

^ We ignore Coma Ber when calculating the average of a and /3 for YCs 
because the limited number of independent data points in its MF results in 
rather large uncertainties on the fitting parameters. 



the value of nic and the clusters' physical properties such as 
total mass, concentration, age or metallicity. 

For example, the dense GCs in our sample span almost two 
decades in metallicity and total mass, but have practically the 
same nic, while all YCs have practically the same metallicity 
but a wide range of m^. values. Therefore, if the environment 
plays a role in determining the shape of the stellar MF this 
effect does not seem to be the dominant one, at least not for 
the clusters in our sample. In fact, taken at face value, the data 
in Table[T]show no sign of such an effect. 

A more interesting possibility to consider in order to ex- 
plain the observed spread in m^ is the role of dynamical evo- 
lution, which could alter the shape of the MF over time be- 
cause of the preferential loss of low-mass stars caused by the 
two-body relaxation process. It is thus worth investigating 
whether our data show any correlation or trend between the 
characteristic mass and the dynamical state of the cluster 

For this reason we have listed in Table[T]the dynamical age 
Tdyn of each cluster, expressed as f /fdis or the cluster age t in 
units of the dissolution time fjis, in turn defined as the time at 
which the cluster has lost 95 % of its original mass. Using N- 
body models, Baumgardt & Makino (2003; hereafter BM03) 
showed that fdis ex t^^^t^~'' where t^h is the half-mass relaxation 
time, fcr the crossing time and x depends on the initial con- 
centration of the cluster Wq. Later, Gieles et al. (2005), using 
the same N-body models, found that fjis can also be expressed 
as a function of the total initial mass of the cluster Mini as 
fdis oc foA/*J,f^ (see also Gieles 2010), where the tidal disrup- 
tion parameter fo depends on the properties of the host galaxy 
and on the cluster's orbit. 

The values of Tdyn listed in Table[T|come from Baumgardt 
et al. (2008) for the specific GCs in our sample and take into 
account the differences that exist between their orbital param- 
eters. As regards the YCs in the sample, they are all located 
in the solar neighbourhood and should therefore experience 
rather similar interactions with the Galaxy. This is consis- 
tent with the results of Lamers et al.'s (2005) study of the 
age distribution of open clusters in the solar neighbourhood 
within 600 pc, based on the catalogue by Kharchenko et al. 
(2005). They conclude that fo = 3.3;!;joMyr, so if we accept 
an uncertainty of '--^l/S on the value of fdis, we can ignore its 
dependence on the specific cluster's orbit. The individual fdis 
values for the YCs in our sample can then be derived follow- 
ing BM03's N-body simulations for clusters with low initial 
concentration Wo = 5 and a distance of /?c = 8.5 kpc from the 
Galactic centre. 

Using BM03's Figure la, we have estimated the initial mass 
Mini of our clusters from their present-day total mass Mtot and 
age f as given in Table[T] While BM03 consider clusters with 
an initial number of stars A^ ranging from 8K to 128K objects, 
for ages younger than ^ 650 Myr, as is the case for all YCs 
in our sample, the ratio Mtot/Mini varies only marginally with 
A^, at most by 10 % and only imperceptibly for f < 300 Myr. 
We have therefore followed the curve for A^ =32K in Fig- 
ure la of BM03 to derive the estimated value of Mini, which 
we have then used to derive the time to dissolution fdis and 
hence Tdyn by means of Figure 3 in that same paper, again 
following the distribution for clusters with Re = 8.5 kpc and 
Wo = 5. For clusters younger than 100 Myr we have assumed 
fdis — 500 Myr but our conclusions would not change if fdis 
were larger or smaller by a factor of two for these clusters. 

It is to be expected that the values of Tdyn derived in this way 
are somewhat uncertain, even though they are as accurate as 
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Fig. 2. — Run of the characteristic mass nic as afunction of the dynamical 
age T;]y[,. Dense and loose GCs are indicated, respectively, with diamonds 
and crosses, with a filled circle providing the average value of each class. 
The core of IC 348 is indicated as an open circle to distinguish it from the 
whole cluster The dashed line is an eye-ball fit to the data. 

one can obtain based on the available information. Lamers 
et al. (2005) point out that the disruption times obtained by 
BM03 are most likely an upper limit because they are calcu- 
lated for tidal disruption in a smooth tidal field without other 
destruction mechanisms. Hereafter we will assume an uncer- 
tainty of about 45 % on fjis, which combined with a typical 
accuracy of ^ 20 % on cluster ages implies an uncertainty of 
about 50 % on Tdy,,. 

The relatively high Tdyn value of the Hyades is consistent 
with the conclusion of Ferryman et al. (1998) that this clus- 
ter has lost about 3/4 of its original mass. Conversely, Prae- 
sepe might not be as dynamically old as the estimated value 
of Tdyn would suggest if it is confirmed that the cluster is in- 
deed the result of a merger of two clusters with different ages 
(see Holland et al. 2000; Franciosini, Randich & Pallavicini 
2003). Finally, a particularly interesting case is that of Taurus. 
Ballesteros-Paredes et al. (2009) conducted a gravitational 
analysis of its orbit and concluded that the Taurus molecular 
cloud must be suffering significant tidal disruption, in spite 
of its young age (1 Myr). We conservatively assigned to it 
■''dyn = 1/2, although the actual value could be larger 

The run of nic as a function of Tdy,, is shown in Figure|2] 
where all clusters are labeled individually, except for the GCs 
(see caption for details). The figure reveals a rather remark- 
able trend of increasing nic with dynamical age over the entire 
range covered by the data. The dashed line is a purely empir- 
ical eye-ball fit to the data with the relationship 
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although an index of 1/2 would still give an acceptable fit. 
Equally remarkable in Figure|2]is the absence of dynamically 
young clusters with large nic but even more striking, consid- 
ering the logarithmic scale of the abscissa, is the lack of dy- 
namically old objects (rdyn > 0.3) with a small nic, even for 
GCs with a total mass of order 10^ M^ . 



There are several effects that can cause the scatter seen in 
Figure|2] around the dashed line, including the uncertainty on 
Tdyn and nic, the role of binaries and incomplete coverage of 
the clusters. As an example of the latter, we show in Figure|2] 
also the data point corresponding to the sole central regions 
of IC 348 (from the MF of Muench et al. 2003). In spite of 
the young age of the cluster, the MF in its core shows strong 
signs of mass segregation, with nic = 0.33 M© as opposed to 
nic = O.12M0 for the cluster as a whole. This may be due 
to primordial mass segregation and/or to very fast dynamical 
evolution (see e.g. Allison et al. 2009), but in either case this 
data point should be excluded from the sample in order to not 
skew the results. In spite of our attempts to select only clus- 
ters with complete radial coverage, however, this information 
might not be completely available or fully correct in the liter- 
ature. This is likely to cause a larger uncertainty on nic than 
what we obtain through the TPL fit to the MF and could ac- 
count for some of the residual scatter that we observe. An ad- 
ditional source of uncertainty on nic comes from the fact that 
all the MFs considered here are system MFs, i.e. we do not 
attempt to correct for the presence of binaries. For dense GCs 
the binary fraction is just a few percent (Davis et al. 2008) 
and is not expected to significantly alter the value of nic, but 
it is very likely higher for loose GCs and YCs (Sollima et al. 
2007; 2010) and its effect should be considered. We will do 
so in a forthcoming paper. Finally, as mentioned above the 
uncertainty on the value of Tdyn can be important and this is 
witnessed by the size of the error bars that we show in Fig- 
ure|2] 

6. DISCUSSION AND CONCLUSIONS 

At least qualitatively, the trend seen in Figure|2]is fully con- 
sistent with the onset of mass segregation via two-body re- 
laxation in a tidal environment, causing the preferential loss 
of low-mass stars from the cluster and hence a drift of the 
characteristic mass towards higher values. Vesperini & Heg- 
gie (1997) showed from N-body models that stellar evapora- 
tion, integrated over the cluster's orbit and further enhanced 
by the presence of the Galactic tidal field, causes a flattening 
of the MF, i.e. a selective depletion at the low-mass end. Their 
models do not reveal a drift of the characteristic mass towards 
higher values, but their IMF in all cases is a pure power-law. 
A drift of nic is visible in the N-body models of Portegies 
Zwart et al. (2001), although it proceeds very slowly and only 
appears when the study of the MF is limited to the regions 
inside the half-mass radius. 

That mass segregation is present in both young and globular 
clusters has long been established observationally. Examples 
include many of the YCs in our sample, such as M 35 (e.g. 
Mathieu 1983; Barrado y Navascues et al. 2001), Pleiades 
(e.g. Converse & Stahler 2008), Hyades (Reid 1993; Perry- 
man et al. 1998), Blanco 1 (Moreaux et al. 2007). As for 
GCs, mass segregation has been observed in all those objects 
where MF measurements exist at various radial distances (see 
e.g. De Marchi et al. 2007 and references therein) and can 
be satisfactorily explained as being the result of the two-body 
relaxation process (Spitzer 1987). Furthermore, De Marchi et 
al. (2007) have also shown that dense and loose GCs have to- 
day systematically different GMFs and this is consistent with 
a much stronger selective loss of low-mass stars in the latter 
Loose GCs have systematically larger nic than denser clusters 
with otherwise similar metallicity or space motion parame- 
ters. This strongly suggests that the value of nic must drift as 
a result of the cluster's dynamical evolution, which is quicker 
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in loose systems, even though the analysis of the N-body mod- 
els (e.g. Portegies Zwart et al. 2001) has so far failed to detect 
this effect. In fact, Kruijssen (2009) has recently developed a 
simple physical model for the evolution of the MF that builds 
on the pioneering work of Henon (1969). His results (see Fig- 
ure 14 in Kruijssen 2009) show a shift in the inflection point 
of the MF with evolutionary time, consistent with a drift in nic 
(D. Kruijssen, priv. comm.). 

A potentially very interesting consequence of the trend 
shown in Figure|2] if dynamical evolution is indeed at its ori- 
gin, is that it would seem plausible that all GCs, having now 
rric — 0.33 M© if they are dense or ^ 0.65 M© if they are 
loose, were born with a stellar MF like the one measured to- 
day in the youngest Galactic clusters, or more precisely with 
nic ~ 0.15 M0. Although not necessarily required by Fig- 
ure|2l this scenario is consistent with it. 

The careful reader will notice that this statement appears at 
odds with the conclusions of a previous work by our team. 
From a study of the GMFs of twelve halo GCs with widely 
different dynamical histories and noting that they show a very 
similar characteristic mass (m^ = 0.33 ±0.03 M© ), Paresce & 
De Marchi (2000) concluded that this value of nic had to be a 
feature of the stellar IMF in GCs. However, it is now apparent 
that this conclusion was the result of a selection effect, due 
to the fact that all twelve GCs in the Paresce & De Marchi 
(2000) sample have high central concentration, since no re- 
liable information on the low-mass end of the MF of loose 
GCs was available at that time. The analysis presented by De 
Marchi et al. (2007) suggests that dense GCs have lost less 
low-mass stars than looser clusters throughout their life, but 
this does not mean that they have lost none. Therefore, one 
can no longer conclude that the m^ ~ 0.33 M© value observed 
in the present-day GMF of dense GCs reflects the character- 
istic mass of their IMF, since it could have been smaller 

The possibility that GCs were born with a much smaller 
nic value than what we measure today in their MF would 
also be consistent with the absence of a turn-over in the 
MF of the Galactic bulge down to the observational limit at 
^ 0.2 M© (Zoccali et al. 2000). Regardless as to whether the 
bulge is primarily the result of fast formation at early epochs 
(e.g. Ballero et al. 2007) or a collection of the stars lost 
from disrupted clusters throughout the life of the Galaxy (e.g. 
Gnedin & Ostriker 1997), its very location at the bottom of 
the Galaxy's potential well makes it very hard for low-mass 
stars to escape from it. Therefore, while the present MF of the 
Galactic disc could represent a complex average over time and 
space, that of the bulge should have not been altered by dy- 



namical evolution and, unlike the case of the comparably old 
GCs, should still reflect the properties of the IMF. Forthcom- 
ing HST observations of the Galactic bulge (see e.g. Brown 
et al. 2009) will provide insights on its MF below 0.15 M©, 
where we expect a turn-over. 

To be sure, N-body models of dense stellar systems able 
to deal with a large number of particles (A^ > 10^) and incor- 
porating a realistic treatment of stellar evolution and of the 
interaction with the Galactic tidal field (see Portegies Zwart 
et al. 2010) would be necessary to simulate the growth of the 
characteristic mass with time in GCs. Our analysis (see Sec- 
tion 5) suggests that young and globular clusters may share 
rather similar values of the a and /? parameters. Thus, if the 
picture sketched above is correct and, by strongly affecting 
the value of nic, dynamical evolution act as the main source of 
variations in the present-day stellar MF of Galactic clusters, 
one might postulate that Population I and II stars should have 
similar IMFs. While not necessarily required by our analysis, 
this hypothesis is fully compatible with our results. Claims 
of the IMF universality that one finds in the literature (e.g. 
Gilmore 2001; Elmegreen et al. 2008) are based on the obser- 
vation that stellar MFs in different environments have similar 
nic values, within a factor of a few. Here we show that it is 
plausible that the original value of nic be actually the same for 
Population I and II stars. 

If this is true, one would have to consider the possibility that 
the characteristic mass in the IMF is not necessarily set by the 
Jeans mass scale in the forming cloud (e.g. Larson 1998; Lar- 
son 2005) and that the physical conditions of the environment 
do not play a significant role in the outcome of the star for- 
mation process, i.e. on the IMF. This, however, would not 
necessarily mean that the star formation process itself is the 
same: recent numerical simulations proposed by Goodwin & 
Kouwenhoven (2009) and Bate (2009) suggest that the IMF 
shape is largely insensitive to parameters such as the binary 
fraction and mass ratio in binaries. If this is the case, the role 
of the IMF as an effective tool to probe the physics of star 
formation may have to be critically reconsidered. 
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hospitality via the visitor programme. SPZ acknowledges the 
support of NWO through grant # 639.073.803. 



REFERENCES 



Allison, R., Goodwin, S., Parker, R., de Grijs, R., Portegies Zwart, S., 

Kouwenhoven, M. 2009, ApJ, 700, L99 
Ballesteros-Paredes, J., Gomez, G., Loinard, L., Torres, R., Pichai'do, B. 

2009, MNRAS, 395, L81 
Ballero, S., Matteucci, E, Origlia, L., Rich, M. 2007, A&A, 467, 123 
BaiTado y Navascues, D., Stauft'er, J., Bouvier, J., Martin, E. 2001, ApJ, 546, 

1006 
BaiTado y Navascues, D., Stauft'er, J., Bouvier, J., Jayawardhana, R., 

Cuillandre, J.-C. 2004, ApJ, 610, 1064 
BaiTado y Navascues, D., Staufter, J., Jayawardhana, R. 2004, ApJ, 614, 386 
Bate, M. 2009, MNRAS, 397, 232 

Baumgardt, H., De Marchi. G., Kroupa, R 2008, ApJ, 685, 247 
Baumgardt, H., Makino, J. 2003, MNRAS, 340, 227 
Bouvier, J. et al. 2008, A&A, 481, 661 
Brown, T. et al. 2009, AJ, 137, 317 
Caballero, J. 2008, A&A, 478, 667 
Chabrier, G. 2003, PASP, 115, 763 



Chappelle, R., Pinfield, D., Steele, 1., Dobbie, P, Magazzvi, A. 2005, 361, 

1323 
Converse, J., Stabler, S. 2008, ApJ, 678, 431 
Corey, K. et al. 2008, AJ, 136, 1778 
Davis, D., Richer, H., Anderson, J., Brewer, J., Hurley, J., Kalirai, J., Rich, 

R., Stetson, R 2008, AJ, 135, 2155 
De Marchi, G., Paresce, E, Portegies Zwart, S. 2005, in The Initial Mass 

Function 50 years later, eds. E. Corbelli, F. Palla, H. Zinnecker 

(Dordrecht: Springer), 77 
De Marchi, G., Paresce, E, Pulone, L. 2000, ApJ, 530, 342 
De Marchi, G., Paresce, E, Pulone, L. 2007, ApJ, 656, L65 
De Marchi, G., Pulone, L., Paresce, F 2006, A&A, 449, 161 
Dolan, C, Mathieu, R. 2002, AJ, 123, 387 
Elmegreen, B. 2010, in lAU Symp. 266, eds. R. de Grijs, J. Lepine 

(Cambridge: CUP), 3 
Elmegreen, B., Klessen, R., Wilson, C. 2008, ApJ, 681, 365 
Forbes, D., Bridges, T. 2010, MNRAS, 404, 1203 



Evolution of the stellar mass function in Galactic clusters 



Franciosini, E., Randich, S., Pallavicini, R. 2003, A&A, 405, 551 
Garmany, D., Conti, R, Chiosi, C. 1982, ApJ, 263, 777 
Gieles, M., Bastian, N., Lamers, H., Mout, J. 2005, A&A, 441, 949 
Gieles, M. 2010, in lAU Symp. 266, eds. R. de Grijs, J. Lepine (Cambridge: 

CUP), 69 
Gilmore, G. 2001, in Starburst Galaxies, eds. L. Tacconi, D. Lutz 

(Heidelberg: Springer), 34 
Gnedin, O., Ostriker, J. 1997, ApJ, 474, 223 
Goodwin, S., Kouwenhoven, M. 2009, MNRAS, 397, L36 
Harris, W. 1996, AJ, 112, 1487 

Hodapp, K., Iserlohe, C, Stecklum, B., Krabbe, A. 2009, ApJ, 701, LlOO 
Holland, K., Jameson, R., Hodgkin, S., Davies, M., Pinfield, D. 2000, 

MNRAS, 319, 956 
Kharchenko, N. Pisnukov, A. E., Roser, S., Schilbach, E., Scholz, R. 2005, 

A&A, 438, 1163 
Kraus, A., Hillenbrand, L. 2007, AJ, 134, 2340 
Kroupa, P. 2002, in Modes of Star Formation and the Origin of Field 

Populations, ASP Conf. Ser 285, eds. E. Grebel, W. Brandner (San 

Francisco: ASP), 86 
Kruijssen, D. 2009, A&A, 507, 1409 
Kruijssen, J., Portegies Zwart, S. 2009, ApJ, 698, L158 
Kucuk, I., Akkaya, I. 2010, RMxAA, 46, 109 
Lada, C, Lada, E. 2003, ARAA, 41, 57 
Lamers, H., Gieles, M., Bastian, N., Baumgardt, H., Kharchenko, N., 

Portegies Zwart, S. 2005, A&A, 441, 117 
Larson, R. 1998, MNRAS, 301, 569 
Larson, R. 2005, MNRAS, 359, 211 
Luhman, K. 2007, ApJS, 173, 104 
Luhman, K., Rieke, G. 1999, ApJ, 525, 440 
Mathieu, R. 1983, ApJ, 267, L97 
Moraux, E., Bouvier, J., Stauffer, J., Barrado y Navascues, D., Cuillandre, 

J.-C. 2007, A&A, 471, 499 



Moraux, E., Bouvier, J., Stauffer, J., Cuillandre, J.-C, 2003, A&A, 400, 891 

Muench, A. et al. 2003, AJ, 125, 2029 

Paresce, R, De Marchi, G. 2000, ApJ, 534, 870 

Pen-yman, M., et al. 1998, A&A, 331, 81 

Piskunov, A., Schilbach, E., Kharchenko, N. V., Roser, S., Scholz, R. 2008, 

A&A, 477, 165 
Portegies Zwart, McMillan, S., Gieles, M. 2010, ARAA, in press 

(astro-ph/1002.1961) 
Portegies Zwart, McMillan, S., Hut, P, Makino, J. 2001, MNRAS, 321, 199 
Reid, N. 1993, MNRAS, 265, 785 
Reid, N., Gizis, J., Hawley, S. 2002, AJ, 124, 2721 
Richer, H. et al. 2004, AJ, 127, 2771 
Richer, H. et al. 2008, AJ, 135, 2141 
Salas, L., Cruz-Gonzalez, I. 2010, RMxAA, 46, 37 
Salpeter, E. 1955, ApJ, 121, 161 
Scalo, J. 1986, Fund. Cos. Phys., 11, 1 

Slesnick, C. Hillenbrand, L., Carpenter, J. 2004, ApJ, 610, 1045 
Sollima, A., Beccari, G., Ferraro, F., Fusi Pecci, R, Sarajedini, A. 2007, 

MNRAS, 380, 781 
Sollima, A., Carballo-Bello, J., Beccari, G., Ferraro, F, Fusi Pecci, F., 

Lanzoni, B. 2010, MNRAS, 401, 577 
Spitzer, L. 1987, Dynamical evolution of globular clusters, (Princeton: PUP) 
Van Buren, D. 1985, ApJ, 294, 567 
Vesperini, E., Heggie, D. 1997, MNRAS, 289, 898 
Vesperini, E., McMillan, S., Portegies Zwart, S. 2009, ApJ, 698, 615 
WiUiams, D., Rieke, G., Stauffer, J. 1995, ApJ, 445, 359 
Zoccali, M., et al. 2000, ApJ, 530, 418 



